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Os números e o computador 


META 


Associar os conceitos de algoritmo, representação dos 
números no computador e implementar cálculos 
usando algoritmos. 


OBJETIVOS 


Identificar os tipos de algoritmos, as propriedades e a 
forma de armazenamento na memória do computador. 
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1.1 Introdução 


Com o aparecimento dos computadores na década de 40, muitos problemas 
foram resolvidos através da aplicação de métodos numéricos, o que antes sem a 
utilização das máquinas eram inviáveis pelo grande esforço de cálculo manual. 

Os homens, através do tempo, preocupam-se com formas de facilitar os 
cálculos, exemplos: o ábaco, inventado pelos babilônios, e os kipus inventado pelos 
incas. 


Kipus 


A tecnologia dos computadores foi avançando cada vez mais, em termos de 
exatidão e tempo de execução das instruções, porém as propriedades da aritmética 
Real não são válidas quando são executadas no computador, pois a memória do 
computador é finita. 

A matemática aborda estes problemas no ramo da Análise Numérica e 
Cálculo Numérico. 

A UFS como outras universidades federais, tiveram computadores desde os 
anos de sua fundação, do tipo IBM, 1130, 360, 3090, e computadores pessoais 
Cobra, Itautec e outros. Só para ter uma idéia o IBM 1130 tinha 8 Kb de memória 
Ram, mas nessa configuração, rodava-se o sistema acadêmico, folha de pagamento 
e vestibular, claro o número de usuários era bem menor. Epocas do cartão perfurado, 
a linguagem utilizada era o FORTRAN Comercial. 


IBM 360 IBM 3090 


IBM 1130 


O estudo da matemática pode ser visto sob dois grandes aspectos: 


Matemática Pura e 
Matemática Aplicada. 


Dentro da matemática Aplicada encontra-se a Matemática Computacional. 
Esta usa como ferramenta o computador e utiliza-se da Teoria da Computação, da 
Teoria da Informação e da Teoria dos Algoritmos. 

A matemática computacional pode ser dividida em três áreas: 


Matemática Simbólica 


Matemática Gráfica 
Matemática Numérica 
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A matemática Simbólica trata dos dados em forma literal, obtendo uma 
solução exata não numérica. É também chamada de matemática não-numérica. Por 
exemplo, provar teoremas utilizando o computador para construir ou verificar 
sequência de inferência lógica que conduzam à demonstração. 

A matemática Gráfica trabalha com dados de forma gráfica e o resultado 
também é um gráfico. As aplicações podem ser divididas em três áreas: 
Processamento de imagens, Reconhecimento de Padrões e computação Gráfica 
Gerativa. 

A matemática Numérica trata da solução de problemas matemáticos através 
do computador e dar como resultado aproximações numéricas. Ela engloba várias 
disciplinas, tais como: Cálculo Numérico, Análise Numérica, Aritmética 
Computacional, Álgebra Numérica, Estatística Numérica, etc. 

O Cálculo Numérico usa métodos construtivos para a solução dos 
problemas, e utiliza só operações aritméticas elementares ( +, -,*,/) e através delas 
são implementadas as demais operações mais complexas. A forma como é 
implementada no computador, o processo de cálculo para a solução do problema 
denomina-se Algoritmo. 

Algoritmo, em geral, é uma sequência finita de passos e operações 
ordenadas que levam à solução de um problema. 

Os algoritmos podem ser numéricos e não numéricos. Os algoritmos 
numéricos são aqueles que utilizam operações aritméticas. 

Exemplos de algoritmo não numérico: 


e Uma receita de bolo 
e Trocar um pneu de um carro 
e Construir uma casa 


Exemplos de algoritmo numérico 


e  Multiplicar duas matrizes Anxp * Bpxm 
e Calcular o sen(x) por uma soma de Taylor 
e Calcular as raízes de um polinômio de grau 2 


Um algoritmo de boa qualidade deve ter as seguintes características: 


Inexistência de erro lógico 

Inexistência de erro operacional 
Quantidade finita de cálculos 

Critério de exatidão 

Independência da máquina 

Os limites do erro devem convergir a zero 
Eficiência 


ESuADa OA e EAD DA 
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1.2 Tipos de Algoritmos 


Algoritmo por computação discreta 
E obtido por uma sequência de computações elementares. 
Exemplo: Algoritmo de Báskara 


Algoritmo por enumeração 
É o tipo de algoritmo que experimenta todas as possíveis respostas em uma 
certa ordem para encontrar a melhor solução do problema. 
Exemplo: Busca do melhor caminho em grafos. 


Algoritmo iterativo 
O algoritmo iterativo ou repetitivo encontra uma série de respostas 
aproximadas que gradualmente vão se aproximando da resposta cor reta até que um 
critério de parada seja atingido por exatidão ou número de repetições. 
Exemplo: Gera em forma repetitiva uma sequência de valores, que 
deverá se aproximar à solução. Terá essa sequência uma propriedade de 
convergência. Algoritmo de Newton para zeros de funções. 


Algoritmo por divisão e conquista 

O problema é dividido em vários subproblemas do mesmo tipo, mas 
menores, que podem ser resolvidos diretamente ou subdivididos novamente, usa-se 
esta técnica, até que todos os subproblemas possam ser resolvidos. 

Exemplo: Dado um intervalo [a,b] onde uma função continua troca de 
sinal, f(a)*f(b) < O. Encontrar o f(c)=0. O algoritmo da bisseção divide o intervalo na 
metade e verifica nas partes onde continua trocando de sinal, descartando a outra 
parte. 


Algoritmo por tentativa e erro 
Este algoritmo procura uma possível solução(tentativa). Caso esta não 
seja (erro), volta à busca segundo novos critérios. E assim por diante. 
Exemplo: Encontrar o menor n tal que (n2 +1)/n! < 10-5, para n Z, 


Algoritmo guloso 
Este algoritmo é usado em problemas de combinatória, onde se busca uma 
solução rápida. Em um processo de escolha sempre é eleito o mais barato. 
Exemplo: Na busca de caminhos em grafos por camadas, em cada 
expansão de um nó escolhe-se o menor se é custo, ou o maior se é lucro. A solução 
não é a melhor, mas ela é sub-ótima. 
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1.3 Solução de um Problema usando o Computador 


Para resolver um problema utilizando o computador devemos seguir pelo 
menos as seguintes etapas: 


1. Selecionar a área onde se encontra o problema no mundo real. 

2. Formalizar o problema, levantando informações relevantes ao 
sistema, com a finalidade de estabelecer um modelo que se 
aproxime ou simule tal problema. 

3. Modelação do problema, neste nível deve ser feita a abstração 
dos dados, identificando-se objetos, operações, e variáveis. 

4. Escolha do algoritmo eficiente, definindo a estrutura lógica do 
algoritmo. 

5. Implementação do algoritmo, neste nível escolhe-se a máquina e 
a linguagem a ser utilizada para a definição física e programação 
do algoritmo. 

6. Validação dos resultados. 


1.4 Representação dos números no computador 


O computador é construído ao redor de uma Unidade Central que 
compreende: 


Memória Central 

Unidade de Cálculo (aritmética e lógica) 
Unidade de Controle 

Unidades de entrada e saída 


A memória central está composta de um conjunto de células elementares 
idênticas (BIT), agrupadas em número de capacidade fixa; cada grupo representa um 
BYTE de oito pontos magnéticos; estes grupos estão numerados de zero an, e o 
número de cada um deles é denominado endereço. 

O armazenamento dos números é feito nestes grupos de células. Numa 
máquina digital, cada célula tem dois possíveis estados, que podem ser 
representados como positivo e negativo, ligado ou desligado e O ou 1. 

Na memória do computador cada número é armazenado em um conjunto 
fixo de bits, sendo o primeiro o bit de sinal. 


DORADADA 


ste lides 


wood Mé-cits 3 uyien) 


Bite Byte 


Existem dois modelos de representação dos números: Ponto fixo e Ponto 
flutuante. 
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O sistema de ponto fixo é representado em duas partes; uma parte inteira e 
outra fracionária em uma certa base numérica. A notação é P(b,n,f) onde b é a base 
utilizada, n o número de dígitos da parte fracionária f. 

Este sistema não funciona para números muito grandes ou muito pequenos. 
Por exemplo, o número de Avogadro 0.60225x1023 moléculas por mol não poderia 
representar-se neste sistema. A implementação desta representação como produto 
de uma fração e potência de 10 está no sistema de ponto flutuante. 


Definição 1. - Um número X que é representado em ponto flutuante tem a forma: 
X=mbe;m (-1,-0.1][0.1,1);eZ 
onde: 


m é a parte fracionária cnamada mantissa 
b é a base numérica utilizada 
e é o expoente ou característica 


Exemplo: 
Se N=234,789, X= 0,234789 x 103 ou X=0,00234789 x 105 


Definição 2. - Um número ponto flutuante está na forma normal (normalizado) se o 
valor da mantissa m pertence ao intervalo (-1,-0.1] [0.1,1). 

Exemplo: 

Se X = 0.0154 x 10-2, a forma normalizada é igual a 0.154 x 10-3 


Definição 3. - Diz-se que um número representado em ponto flutuante está na 
“Forma Standard t dígitos na mantissa” se ele está normalizado e esta mantissa tem 
exatamente t dígitos. Se a mantissa tiver mais de t dígitos então arredondar o dígito 
t+1 assim: 

Se o (t+1)-ésimo dígito for igual ou maior que 5 então o t-ésimo dígito é 
incrementado em 1. 

Em outro caso, os t primeiros dígitos são considerados como mantissa. 

O sistema de ponto flutuante é representado por F(b,t,e1,e2) onde: 


b é a base 

to número de dígitos na mantissa 
e1 é o menor expoente 

e2 é o maior expoente 


Exemplos: 

1. Seja t=8, N=0,8934572834 o número está normalizado t dígitos na 
mantissa como X=0,89345728 x 100 

2. Seja t=5, N=3,14159 então X=0,31416 x 10 
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Como se pôde observar nestes dois casos os números originais não 
puderam estar representados completamente para esses computadores hipotéticos, 
por tanto os valores armazenados são aproximados. 


1.5 Resumo 


Nesta aula, você verificou que o armazenamento dos números nem sempre é exata 
quando eles são transformados para uma base diferente da decimal. Isto acarreta 
uma aproximação, por tanto há existência de erro, assunto que veremos na próxima 
aula. 


1.6 Atividades 


1. Seja o número N = 56783945783245 e um computador com t=8, qual é a 
representação em ponto flutuante, precisão simples, e em precisão dupla”? 


2. Verifique se as duas expressões a seguir podem ser usadas para calcular a 
abscissa da interseção da reta,que passa pelos pontos (x0,y0) e (x1,y1), com o 
eixo x. x=(x0y1 — x1y0)/(y1-y0) e x=x0- [(x1-x0)y0]/(y1-y0) 


3. Usar os pontos (1.31, 3.24) e (1.93,4.76) e t=3 dígitos, calcule o x usando as 
fórmulas do exercício 2. Comente. 


4. Cada computador tem o “t” número de dígitos que trabalha, o algoritmo que segue 
calcula este número (resultado em j) 
Pt. e=1 
P2. j=i 
P3. Enquanto 1+e > 1 
j=i+1 
e=e/2 
Fim enquanto 
P4. Mostrar j 


A implementação na linguagem do Scilab 


e=1; 

j=1; 

while 1+e > 1 
j=i+1; 
e=e/2; 

end 

j 
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1.7 Comentário das atividades 
Os exercícios 1 a 3 é para fixar as definições de ponto flutuante. 


O exercício 4 é interessante para descobrir o número t e também um ótimo exercício 
para iniciar na programação no software SciLlab. 


1.8 Referências 


CUNHA, Cristina. Métodos Numéricos. 2º Ed. Campinas SP: Editora da UNICAMP, 
2008. ISBN: 85-268-0636-X, CDD — 620.00151 
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Erros 


META 


Conceituar o erro, as fontes e formas de expressar 
estes erros, propagação dos erros em operações 
aritméticas fórmula geral e problema inverso. 


OBJETIVOS 


Resolver problemas práticos de erros em funções de n 
variáveis e calcular a cota para o erro em processos 
infinitos. 


2.1 Erros 
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Na aula anterior vimos que nem sempre a aritmética computacional coincide 
com a aritmética real. 

Por exemplo, o número Tr, número irracional com infinitos dígitos, não é 
possível ser armazenado na memória do computador por ela ter tamanho finito e fixo. 
Ao armazenar este na forma standard, ocorre uma aproximação do valor exato. Logo 
existe um erro. Neste caso, de arredondamento. 


2.2 Tipos de Erros 
Pela fonte onde são produzidos estes erros podemos classificá-los como: 


Erros de modelação 

Erros inerentes aos dados de entrada 
Erros de arredondamento 

Erros de truncamento 


Os erros de modelação são provenientes da simplificação das situações 
reais que se faz através do modelo, ignorando-se certos aspectos do mundo real. 

Os erros inerentes são os erros cometidos nos valores dos dados, causados 
pela inexatidão das medidas tais como distância, tempo e temperatura, e que são 
medidos por instrumentos limitados, por enganos pessoais ou pela natureza. 

Os erros de arredondamento são o resultado da representação de um 
número numa máquina. 

Os erros de truncamento são erros cometidos pela aproximação de um 
cálculo infinito por outro finito. 


Definição 4. - O erro absoluto é definido como a diferença do valor exato e valor 
aproximado. 


e=Ve-Va 
onde: 
e é o erro absoluto 
Ve é o valor exato 
Va é o valor aproximado 
O módulo do erro absoluto (erro absoluto máximo) é o valor absoluto de: 
jel=|Ve-Va| 
O erro relativo é representado como o erro absoluto dividido pelo valor exato 


dO=€e/Ve = e/Va 


Obs. Usa-se o valor aproximado Va se não se conhece o valor exato Ve 
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O erro percentual é representado como: 


P = 1006 


2.3 Propagação do erro nas operações aritméticas 


2.3.1 Erro na soma e diferença 


P1. Seja Vas = Va; + Vas 

P2. Seja Ves = Ve; + Ve> 

PS. Ej= Ve; E Vai 5 i= 1,2,3 

P4. Vas + €3 = (Var + €1) + (Vaz + €2) 
Pb.es=€ +e> 

P6. le] = |e1 + €4] 

P7. les] < e1| + |ea| 


2.3.2 Erro no produto 


P1. Seja Vas = Vas. Vaz 

P2. Seja Ves = Ve; . Ve 

P3. Ej= Ve; = Vai j i= 1,2,3 

P4. Vas + €g= (Var + €1) . (Vaz + €2) 

P5. Vas+es=€1.e0+€4. Vaz+ Va. eo + Vas Vas 
P6. Es = €1 «co + Vas. €4 + Va. €o 

P7. €3/ Vas = (Vas. € + Va; (Es +E. €2) / Vas 

Ps. Os = 01 .0o+ 01 + O2 

P9. |ôs| = [of E Oo + Õ1 + d2| 

P10. |ôs| < |04] . 02] + |04] + |d2] 


2.3.3 Erro na divisão 


P1 Seja Vas E Va; / Vas 
P2.Seja Ves = Ve; / Ves 
P3E= Ve; = Vai ; |= 1,2,3 
P4.Vas + €3 = (Va; + €1) / (Vas + €2) 
P5.Vas+es=(Val+e).(1/Vaz).(1+€e2/Vao)” 
P6.Vas + €3= (Vas/Vaz+tes/Vas). (1-7...) 
P7.e3= (Va2.€1- Val .€2)/ (Vaz) 
P8.€e3/Va3=(1/Va3). (Va2.€1-Val.e2)/ (Vaz) 
P9.03=51-52 
P10./53 | = |51 - 52] 
P11./53] < |51] + |52] 

Exemplo 1: 
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(Hipótese) 
(Hipótese) 
(def. de erro absoluto) 


(definição de módulo 
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Sejam os números irracionais m = 3.1415926535897931159980...e 
V2 = 1.4142135623730951454746 ... (valores dados pelo Scilab , %pi e sart(2) com 
format(25)), que devem ser armazenados em uma máquina de t=8 dígitos na 
mantissa. Qual é o erro absoluto para t+ v2 e para mr + 2? 


Os valores de 17 = 0.31415927 + 10! e V2 = 0.14142136, estão na forma 
standard normalizada. 

Observe que as mantissas foram arredondadas. Pela definição de erro 
absoluto: 


Er = 3.1415926535897931159980... — 3.1415927 = —0.0000000464102067887495.... 
E = 14142135623730951451746... — 14142136 = —0.0000000376269049251476 ... 


Ema = En teya= —0.0000000464102067887495... — 0.0000000376269049251476... 
= —0.0000000087833018636019... 


Emi) = Tre + VZ*E+ Eres 
= 3.1415927 « (—-0.0000000376269049251476) + 1.4142136 
* (—0.0000000464102067887495) + (—0.0000000376269049251476) 
* (—0.0000000464102067887495) = 0.0000927299886941268531 


Exemplo 2: 
Qual o erro máximo para um número x com t dígitos na mantissa, se ele é 
arredondado? 


O dígito t é acrescentado uma unidade se o dígito t+f > 5, em outro caso 
não se modifica. Então os valores do dígito t e dígito t+1 poderão sertOt1t2t3out 
4 ou (t-1)5, (t+1)6, (t+1)7, (t+1)8 ou (t+1)9, fazendo parte do valor exato, e os 
respectivos erros £0,€1,€2,€3,€4,€5,€6,£7,£8,£9 Serão menores que 0.5x107. 


2.4 Fórmula geral para Erros 


Seja uma função diferenciável y = f(x). 
Uma variação de x em um (x) faz que mediante a função f haja uma variação 
de y, assim: 


A(y) = T [x + 09 ]- f69), 


considerando A(x) um valor muito pequeno, ou seja, uma quantidade que 
represente um erro em x. Utilizando f teremos uma variação de y, chamada de A(y), e 
que representará o erro em y ou da função f. 

Se em A(y) = f[x + A(x) ] - f(x) multiplicando e dividindo o segundo membro 
por A(x) temos: 
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A(y) = f(x) Ax, considerando Ax muito pequeno 
ley |= [FCO] lex 


Em geral seja: 
U=f(Xy,X2, X3,..., Xn) 
seja: 
Ex parai=1,2,3,..,n 
os erros absolutos dos argumentos de f, então o erro absoluto de u é: 
Eu=(0f/ 0x1) ex + (Of / Oxo) Exo +... + (Of / Oxn) Exn 
leul = I(0f/ 0x1) ex + (Of/ xo) Exo +... + (0f/ 0Xn) Exa] 
leu < I(0f/ 9x1)] |ex1] + I|(9f/ 0x5)] |ex2] +... + |(0f/ 0xn)] lex] 
du <I(Of/ 9x/U) ex] + I(Of/ Oxo/U) Exo] +... + |(0f/ Oxn/U) Exnl 
Exemplo: 
Achar o máximo erro absoluto e relativo do volume de uma esfera se o 
diâmetro D = (3,7 + 0,05) cm, 11 = 3,14. 
V=1/6M Dº=vV(n,D) (função em duas variáveis) 
Observe que é considerada uma variável porque tem erro. 
€p < 0,05 En < 0,0016 


oV/0M = 1/6 D*= 1/6 (3,7)º = 8,4421666 
oV/9D=1/2ND2= (1/2)3,14 (3,7)? = 21,4933 


ev =(9V/0Men+ (0V/ 9D)ep = 

< 8,4421666 (0,0016) + 21,4933 (0,05) = 

< 0,01350746656 + 1,074665 = 1,08817246656 cm? 
dv < 1,08817246656 / 26,508403 = 0,041050 ou 

dy < (1/6 Den +1/2M Dº ep)/(1/6 MD?) 


Ov<On+30p 
dv < (0,0016/3,14) + 3(0,05/3,7) = 0,0005095 + 0,0405405 = 0,04105 


2.5 Problema inverso do Cálculo de Erros 


O problema inverso do cálculo de erros consiste em encontrar os erros dos 
argumentos de uma função dado o erro da função. 


Seja a função: 
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U=T(X1, X2, X3, X4, ..., Xn) 


dado o erro em u determinar os erros para X1, X2, X3, X4, ...; Xn- 

Este problema não tem solução analítica exata, já que temos n incógnitas 
para uma única equação. 

Para poder dar uma solução a este problema há que se considerar 
restrições a fim de reduzir o problema a uma equação e uma incógnita. 

Uma alternativa para esta redução é considerar “O princípio de Efeitos 
iguais”, hipótese que usa-se em estatística. Este princípio supõe que as leis físicas 
atuam da mesma maneira para uma ação produzindo efeitos iguais. 

Extrapolando esta idéia para o problema inverso do cálculo de erros surgem 
as seguintes hipóteses: 


|. Os erros absolutos são iguais para x1, x2, ..., Xn. 


Il. Os erros relativos são iguais para x1, X2, ..., Xn. 
III. A função u contribui no erro total. 


2.5.1 Hipótese | 


P1. Exi = Ex2=Ex3=...=yn= ki (Hipótese) 
P2. ey é conhecido (Hipótese) 
PB. p= (01/03) Easl= 12, aa (def. de £) 
P4. eu=ky D(0f/ Oxi) (P3, P1) 
P5. ki =eu/X(0f/ Oxi) (P4) 


2.5.2 Hipótese II 


P1.0m =0w=...=Om=kp (Hipótese) 
P2. ey é conhecido (Hipótese) 
P3. Eu= L(0F / 08) Ea sl= 1,2, o; (def. de £) 
P4. cu=D(0f/0Oxi) ex xi/xi,i=1,2,...,n (P3) 
P5.eu=5(0f/ôx)xida,i=1,2,...,n (P4) 
P6. cu=k23(0f/0Oxi)xi,i=1,2,...,n (P5, P1) 
Pr = eu 2(0f/ 0) = 2, sas (P6) 


2.5.3 Hipótese III 


P1.(0f/0x1) x1 = (0f/ 0x2) x2 =... = (0f/ dxn) xn =k3 (Hipótese) 
P2. ey é conhecido (Hipótese) 
P3. cu=5(0f/0xi) Ex, [=1,2, ss; À (def. de €) 
P4. cyu=nkg (POP) 
P5.kg=€/n (P4) 


2.6 Erros de Truncamento 
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Erro cometido ao aproximar um cálculo infinito por outro finito. 
Exemplo 1: 


Xo; X1, Xo; X3; mona ig Xn; Xn+1; ne 
X1 = F(xo) 
X2 = F(x1) 


Xn41 = F(Xn) 
Ka 


Seja x uma solução aproximada da sequência ( xn ) que converge a x* no 
limite. 


x > valor aproximado 
x* > valor exato 


RW ehmsa dl 
Xo, X1, X2; «1, X, =. SP X 
c=Vo-Va=x*-% 


* 


onde € é o erro de truncamento. 


Exemplo 2 
Seja uma função f(x), n vezes contínua e diferenciável, num intervalo [a,b]. 


fQ) = Lã Fe) (x— a) Série de Taylor 
se a = 0, Série de Mac-Laurin 


Cálculo da função sen x, isto é, expressar como uma série de potências 


Solução: 
f(x) = sen x f(x) = -cos x 
f(x) = cos x f(x) = sen x 
f(x) = -sen x f””(x) = cos x 


função trigonométrica cíclica 


Paraa=0: 
f(0) = O f(x) = -1 
f(O) = 1 f”(x)=0 
f(0)=0 Tg) = 


0 1 0 1 0 
= = 0 pái J 1 ne J 2 o nr 3 a J 4 Rm 
senx = qi (% 0)º + T (x — 0) + te 0) tal 0) +a (x — 0)* + 
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E 
senx =x 31 + E 
x2i+1 j . 
Ssenx= Le Emo 1)', x radianos 
= do eram 1)! = Ve 
k x2i+1 = io 
si ri ( Di= Ve 
x2i+1 i 
e=Vo-Va= ba =k+1 (+41)! ora Ú 


Estima o erro — Encontrar uma cota superior para o erro. 


Teorema: “Seja a série So(—1)' , tal que JUo| > |U+| > |uz] > |us] >... > |un] > [Uno] >... 
alternada e convergente, então 5 29(—1)' < |uks1| 


x2i+ 


Para a sériesemx:u= - 0< x<TI/4 
(2i+1)! 
g2i+1 x2i+3 x2i+1y2 
—————— | > — | 
(Qi + 1)! 2i+3! Qi+ D!Qi+BMQIi+3) 
x2i+1 x? g2i+1 
(Qi+ D!| |Qi+2D(i+ 3| S 
Para a função sem x: 
x2k+3 
(Qk + 3)! 


Exemplo: 

Determinar o erro de truncamento para o cálculo de sen 30º pela fórmula de 
Taylor. 

Solução: 


“Cate ator 


mm (63, (m/6(m/6 


sen 30º = sen — 


6 6 6 6.20 
x” o (m/6)?” 

SESI Rê 
So = A: 


linoOufs.br 


CÁLCULO NUMÉRICO Aula 2 


g E x2 i 
E 0324 ) 


2 


E ne dd so 1 

ao as. 1756 ) 

E x x? aê x? : 
3 51'6.7 2674 ) 


2 


XxX 
dn = Sr tapeag + DO 


Cálculo do e: 
f(x) = e 


Fazerx=15f(1) = 
f(x) = eS f(x) = Aa 


Parax=1>e=1+% reto to ET + “=>1+1+0,5+ 0,1666 


e = 2,6666 
er=h-U=D 5 
i=k+1 e 
xi at+1 at+2 qt+3 
AE TRID TD RD 
at+1 x x? 
ql PA 
| “RR2 (EDUC) 


x? 


“RIDE FIRTS * 
1 
CERICREINCE CE 
1 


e CULDEADRA 
t+ Ermltr x2 x3 + 


“tRrDi] TE+2 NICE Curas” 

t+ x aqi 1 

E, O | — 

FD ES) x qa USA 
i=0 i=0 

Se a = -<1 
k+2 

x < k+2 
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Se k=0, x<2 
qh+1 1 
Bud Sa 
kK+DUq4. A 
( q rr 
sex= 1, 
Sd k+2 000 k+2 
“CGE D EFZA EAD 
je 5 5 
Er < = 


414 244 96 


Exemplo: 
Quantos termos da série de Taylor são necessários para que o número “e” 
tenha um erro menor que 0,01. 


Solução: 
ET < 0,01 
< ad < 0,01 
ES CER E) A 

para k = 2: Em 

313 18 
para k = 3: Em E 

414 96 
paraks=s4— == 

515 100 


O número de termos é 6. 


2.7 Atividades 


1.Se 1000 aproxima x com erro menor que €. Mostre que 1/1000 aproxima 1/x com 
erro absoluto menor que |£|/| (1000 + c )*| 


2. À altura H e raio R da base de um cilindro são medidos com aproximação de 
0.5%. Qual é o máximo erro absoluto e relativo ao calcular o volume. 7 = 3.14 


3. Um cilindro de alumínio com diâmetro da base d = 2em + 0.01cm altura h = 11cm 
+ 0.02cm e peso p = 93.4gf + 0.001 gf. Determinar o erro relativo do peso 
específico pe = p/v 


4. Determine a série de Taylor para a função f(x) = sen x, e a cota do erro de 
truncamento. Quantos termos da série são necessários para cometer um erro 
menor que 0.01. x = 0.8 


5. Um cilindro de alumínio com diâmetro da base d = 2em + 0.01cm altura h = 1tem 
+ 0.02cm e peso p = 93.4gf + 0.001 gf. Determinar o erro relativo do peso 
específico pe = p/v. 
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Zeros de Funções 


META 


Resolver o problema: dada a função f(x), contínua em 
um intervalo I=[a,b], encontrar um x* tal que f(x*)=0. 


OBJETIVOS 


Estudar diferentes algoritmos, encontrar soluções e 
verificar qual é o mais eficiente e em que condições. 
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3.1 Zeros de Funções 


Um problema comum em engenharia, ou em geral nas diversas áreas das 
ciências exatas, é determinar soluções em equações não lineares. Equações estas 
que envolvem funções transcendentes. 


No problema seguinte: 

Um cabo telefônico suspenso entre dois postes tem um peso de q 
quilogramas-força por metro. A tensão T no meio do cabo é obtida pela resolução da 
equação (2T/qy) senh(pL/2T) = S, onde S é o cumprimento do fio, L é a distância entre 
os postes e q o peso específico do fio por metro linear. 


Os métodos numéricos, expostos nesta aula, são utilizados para encontrar 
soluções aproximadas para equações deste tipo de problema. 

Seja a função f(x) = 0. 

Um zero da função é um x* € R tal que f(x”) = 0. 

Para um intervalo | = [a,b], se f(a).f(b) < 0 c /f(c) = 0. 

f deve ser contínua em | = [a,b]. 


fx) 


Solucão x' 


3.2 Método da Bisseção 
O método é baseado no teorema do valor intermediário que diz: 


“Se uma função é contínua no intervalo I=[a,b] e satisfaz a condição 
f(a)“f(b)<O, então existe pelo menos um x*Efa,b] tal que f(x*)=0.” 


O método também denomina-se de pesquisa binária, e consiste em dividir o 
intervalo inicial em subintervalos que contenha o zero procurado. Divide-se o intervalo 
inicial e se descarta o intervalo onde a função não troca de sinal, prosseguindo com o 
intervalo que satisfaz a condição de troca de sinal. 


3.3 Algoritmo 


P1. Dada a função contínua f(x) no intervalo I=[a,b], definir uma tolerância € e verificar 
se f(a)*f(b) < O. 

P2. c=(a+b)/2 

P3. Se |f(c)| < £ então pare. Solução aproximada c. 
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P4. Se f(a)“f(c) < e então 
b=c 
se não 
a=c 

P5. Volta ao P2. 


Exemplo 
Encontrar uma solução para a equação f(x)=2*-4*x, no intervalo [0,1] usando 
o programa a seguir. 


3.4 Programa no Scilab 


function y=f(x) 
y=2"x-4*x 
endfunction 
format(15) 
a=0; 
=: 
[el=input("Digite a tolerancia Ex.0.0001: "); 
c=(a+b)/2; 
fc=abs(f(c)); 
while fc > e 
if f(a)*f(c) < O then 
D=e; 
else 
a=c; 
end 
c=(a+b)/2 
fc=abs(f(c)); 
end 


3.5 Métodos Iterativos para a Solução do Problema 


Um método iterativo é um método repetitivo, que gera normalmente uma 
sequência de valores. 

Para que a sequência de valores seja gerada deve existir uma fórmula 
recursiva. 

Se a sequência gerada nos leva à solução, então é uma sequência 
convergente. 

O método iterativo possui uma fórmula recursiva, uma regra de parada e 
condições de convergência para a sequência gerada. 

A regra de parada é dada pelo valor de £, chamado também de zero 
numérico. E um valor pequeno, nomeado de tolerância, que indica o grau de exatidão 
da solução. E definido por fora. 
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3.6 Método de Iteração Simples ou iteração linear 


Um número p é um ponto fixo se, para uma função dada g(x), g(p)=p 
Exemplo: 


g(x)=xº-2x+2 


os pontos - 2 e 1 são pontos fixos porque 


Dado um problema f(x)=0 , pode-se definir uma função g(x) de tal forma que 
9(x)=x-f(X) e f0g)=x — g(x) 


Se g(p)=p então g(p) - p = 0 e f(p)=0 e p será um zero. 


F(x)=x)-2x+2 


3.7 Algoritmo 


PO. Transformar a função f(x) = 0, tal que g(x) —x = f(x) =0 
P1. Para j=0 

Escolher um valor inicial qualquer x; em [a,b] 
P2.j-j+1 

x; — 9(Xm) 
P3. (Regra de Parada) 

Se |x;— x; | < € então x; é solução aproximada. 

Se não Voltar a P2 


A sequência gerada é Xo,X1,X2,.............. DN uensuanáii 


linoOufs.br 


CÁLCULO NUMÉRICO Aula 3 


Onde x1=9(X0), x2=9(X1), x3=9(X2)........... Xn=9(Xn1) 
Exemplo 
fo)=x*-2x-1 I=[1,2] 
f(1)=1-2-1=-2 
f(2)=8-41=83 
f(1).f(2) < 0 e=0,01 


PO: g(x) / g(x) - x = f(x) 
x*-2x-1=>x=Y2x+1 


Pt:j—O,xo€I=[1,2] 
Xo=1 


P2: x1 = g(Xo) = 9(1) = Y3 
x1 = 1,44334957 


P3: |xy - xo] = |1,44224957 - 1|> 0,01 


P2: xo = g(x1) = 9(1,44224957) = /(1,44224957).2 + 1 
xo = /3,8845 = 1,571973 


P3: |xo - x1] = [1,5719783 - 1,442251] > 0,01 


P2: xa = g(xX2) = 9(1,571973) = 4/2(1,571973) + 1 


x3 = 1,0622 


P3: |xa - xo] = |1,60622 — 1,571973]| = 0,034 > 0,01 


P2: xa = g(xa) = F(1,60622) = 2/2(1,60622) + 1 


x4,= 1,6150 
P3: |1,6150 — 1,6002] = 0,0088 < 0,01 


x4 = 1,615 Solução aproximada 


3.8 Condições de Convergência 
Dvxe lo g(x)e l; 


II) F(x) deve ser contínua no intervalo |; 
III) O valor absoluto da derivada da função 
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Ig) <1 VxEl. 


A sequência ( xn ) gerada pelo algoritmo de iteração simples convergirá à 
solução x*, isto é, lim, so x, = x“ OU (Xn) =X. 

Exemplo: 

Verificar as condições de convergência para g(x) = V2x + 1. 


)vx€e [1,2], então g(x) € [1,2] 
Parax=1:9g(1) = 1,4422 € [1,2] 
Parax=2:9(2)= 1,71 € [1,2] 
Porque a função é crescente no intervalo então g(x) € [1,2] 


2) ias PR g(x) = 9(Xo) 
Não há pontos de descontinuidade no intervalo 


3) g(x) = (2x+1)'8 


=2 Z =2 2 
g'(x) = 3 (Ux + 1)3 (2x + 1)' = 90x + 1)3 = 


3:/(2x + 1)? 


O is RR 
x= 1: g'(1) = 5 = á(z) <! 
ralije- =|. 
X=2: g'(2) = 3425; — (om) <1 
Porque a função derivada g (x) é decrescente no intervalo 


Observação: O intervalo poderá ser relaxado para satisfazer as condições 


3.9 Considerações das condições de Convergência 


) vx€ lo g(x) € I=[a,b] 
f(x) = 0 > 9(x) — x = f(x) 


II) g(x) é contínua 
Condição |) e Il) garantem a existência de solução 


)asg(a) (1a condição) 
2) bz g(b) (1a condição) 
3)0<g(a)-a=f(a) (1, def. de g(x)) 
4)02g(b)-b = f(b) (2, def. de g(x)) 
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5)fla)>0"fb)<0 (3, 4) 
6)3 ce [a,b]/f(c)=0 (5, def. de zero de função) 


ID|g(t)|V<1,xéEl 
Também conhecida como condição de Lipchitz 


Ig) |<L,L<1,vVxE | 
é equivalente a dizer que para dois pontos q.q. | vale o seguinte: 


Xw X1 € lentão |g(xk) - g(X«1)| < L |x«- X«i] 


Digt)|<sL,L<1,vx€El=Ta,b] (Hipótese) 

2) Seja xe x«: pontos pertencentes ao intervalo | (def. x€ 1) 

3) Existem g(x,) e g(x«1) € | (def. g(x)) 
4) Pelo Teorema do Valor Intermediário: 


3 c [x Xi] / g'(c) = “C-SGk+) 


XkT Xx+1 
5) Ig'(o)l = qeu SC) (4,def. de valor absoluto) 
TXR+1 
6) Como ce |, entãog'(c)<L (5,1) 
7) a <Lou lg(xi) — (Xe )] < Llxp — xp] 


As Condições Ill) e Il) garante a unicidade da solução. 


Ngt<L,L<1,vxE | 


A |g'(xk) - 9 (Xk+1)| < L.| Xk - Xk41 | para 9.0. Xk, Xk+1 €| (Hipótese) 
2) Suponhamos que existem duas soluções S1, s2, diferentes. 
3) s1=9(sS1) ” s>=g(s>) (def. de Solução) 
4)s1,8s2€ |, então |g(s1) - g(s2)| <L |s4 - sol(def. eg. à cond. de Lipchitz) 
5) |ss-s2]<L |sj-sa| (3, 4) 
6) Como s; é diferente de s> — s; - s2% 0, logo 1 <L | contradição com 


a hipótese, constante de Lipchitz menor que 1 
17) > S4=So 


3.10 Interpretação geométrica a condição III 
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As três condições garantem convergência da sequência gerada pelo algoritmo 


Provar que: 
limx, = x* 
n5s 
ou 
limlx, — x*| = O 
nos 


1) |xn - x*| 

2) |xn - x*| = |g(Xn-1) - g(x”)| (def. algoritmo e x* solução) 

3) |9( Xn-1) - 9(X9]| < L.|xn4 - x*] (condição 3 equiv.) 

4) |xn4 - x*| = |g(xn-2) - 9(x)] 

5) |xn - X*| < L.|xn4- x*| S L.L.|xn.2- x*] (2,4) 
|Xn - X*| < Ligas x*| 

6) IXn.z- X| = |9(Xn3) - 9(X)] S L.|Xnes - x] 

7) |xn - X*| < L2.|xn-g - x*| 

8) |xn- x*| < L”.|xo - x*| 

9) limns. |xn - x* < |limn L”.|xo - x*| 

10) lim no L”.| x0 - x* | = [xo - x*|. lim noso L” 

AD limnso Lº=0 porqueL<1 

12) limnso |xn-xX*|=0 


3.11 Erro de truncamento para n passos 


Valor exato: x* 
Valor aproximado: xn, n qualquer. 


Xo; X1; Xo, X3; ...3 Xn-1; Xn 


e=Ve-Va=x*-xn 


* 


D|x*-xn)=limmss |Xm-Xn|;Mm>n 

2) |Xm - Xn) = |Xn - Xml] = |Xn = Xn41 + Xn41- Xn+2 + Xn+20 +.. = Xm1- Xm] 
3) IXn = Xn+1] + |Xn+1 E Xn+2] + |Xn+2- Xn+3] +. + |Xm-1 E Xm| 

4) [en - Xnnl = |9(Xn-1) - 9(xm] SL |Xn1- Xn] 

5) |Xn+1 a Xn+2] = |g(Xn) a 9(Xn+1)| SL |Xn - Xn41l & [EE |Xn- ! - Xn] 

6) |Xn+2- Xn+3] = = |9(Xn+1) - 9(Xns2)| SL | nat a Xn+2 | <|º |Xn- à Xnl 
7) Das Ka] E | ad] La E qa A 
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8) Mm |Xn = Xml < lim m>os |Xn-1 e Xml É ia L 
9) |Xn = x*| 5 |x"!- x] Zy dº 
10) xa -x*]s |"! x] LEG! = | XL 


11) Ixn - x*| < | CuttniCota do erro 


3.12 Atividades 


1. Encontre uma aproximação para 25'º com precisão de 10 usando o algoritmo da 
bissecção 


2. Achar uma função de iteração para encontrar um zero diferente de x = 4 de 2X = 
4x 


3. Dar uma cota do erro de truncamento ao usar n iterações no método de iteração 
simples. A cota deve estar en função dos valores | xo - x1| e a constante de 


Lipchitz L 


4. Demonsire que xn,1 = Xn (2 - Kxn) converge a 1/K quando n tende ao infinito 


3.13 Referências 


CUNHA, Cristina. Métodos Numéricos. 2º Ed. Campinas SP: Editora da UNICAMP, 
2008. ISBN: 85-268-0636-X , CDD — 620.00151 


BURDEN, L. Richard, J. Douglas Faires Análise Numérica SP: Editora Pioneira 
Thomson Learning, 20083. ISBN 85-221-0297-X CDD - 515 
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Zeros de Funções 
(continuação) 


META 


Resolver o problema: dada a função f(x), contínua em 
um intervalo I=[a,b], encontrar um x* tal que f(x*)=0. 
Usando o Método de Newton. 


OBJETIVOS 


Estudar diferentes casos do especiais do Método de 
Newton. 
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4.1 O Método de Newton 


O método de Newton pode ser visto com um caso particular do método de 
iteração linear, onde a função g(x) é construída para ser uma função de iteração. Isto 
é, que satisfaça as três condições de convergência. 

Seja a função f(x) uma função contínua em | = [a,b], e f(a).f(b) < O, com f(x) £ 0 
no intervalo 1. 

Construção da função g(x) de tal forma que x — g(x) = f(x) 


) 
2) PTE = O porser f(x) £ 0 
siri JRR) 
4)x—x Sa 0 
5)x=x- PTE = g(x), onde g(x) é igual ao segundo termo 
6) g(x) satisfaz as condições de convergência do método de iteração 


linear. 


4.2 Interpretação geométrica 


-7 8 ad id 
-2 
-3 
f 
80= LD =P) 
Oo) 
f(x) Po? 
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co TOO) 

O FG) 
o FXn) 
Kn+1 no Kn — f'(x ) 
n 


Exemplos: 
1º) Raiz quadrada de um N>0 


Seja f(x) =x? - N 
x2—N 1[x2+N 
2X. 2 


X =X — 
n+1 n Xn 


Exemplo: 
2 =? 
o — 1 


aoroo= + t]= 302? 


1,17 
ecrio= eB)=iEed =5[e | = 1416666 
x, = F(x,) = 


2º) Raiz K-ésima de um N>0 


f(x) = xe — N 
xk-N 1 N 
foda Aa gd E (e a 


4.3 Algoritmo 


PO. Dada a função f(x) contínua em I=[a,b], f'(x) O em [a,b], f(a).f(b) < O, € (tolerância) 

P1. Escolher um xo inicial em [a,b],j) — O 

P2. j <— + » X=X1— Tx )/f (x) 

P3. (Regra de parada) Se |x-x;1] < £ ou If(x;)| < £, então pare. Solução aproximada x; 
Se não volta ao passo P2. 


1.4 Programa no SciLlab 


deff(Iy]=9(x)', y=2"x-4"x' 
deff([z]=dg(x)','z=log(2)*21x-4') 
x0=0; 

x1=0.5 

format(20) 

eps=0.00001 


linoOufs.br 


CÁLCULO NUMÉRICO Aula 4 


while abs(x0-x1) > eps 
x0=x1; 
x1=x0-(g(x0)Vdg(x0)) 
end 

XI 


1.5 Casos Especiais 
Caso 1. 


Não se pode garantir que a função f(x) seja diferente de zero em todo o 
intervalo i=[a,b]. 

Seja f (xo) £ O, xo E [a,b] 

Considerar f (xo) = M constante para todo o cálculo da sequência. 


= fXn) 
Xn41 = En — f(x) 
f(x0) 0 vn 
(xo) = M £ O 
FXn) 


Kn41 = Xn — M 


Caso 2. 


A derivada da função é complexa. 
Aproximamos a derivada pelo quociente do limite. 


f(xn) = fQn-a) 


F Ca) = ERR Xn — Xn4 
l o f(xn) = f(xn-1) 
f (Xn) ad era 
E RO = TOD 
Xn — REn-1 


Fórmula também conhecida como método da secante. 


Caso 3 
Método de Newton aplicado a polinômios 
Seja o polinômio Ph(x) = ao + aj.X + a2.X2 +... + anXº 


= Pr(xi) 
Xk+1 Es Xk — Pi(xy) 
n 
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Para cada iteração, é necessário calcular o valor numérico do Pn(x) e Ph (xy). 


Exemplo: 

Ps(x) = do + d1.X+ ao.x? + as.xº + ae + as.xº 
Ps(r) = ao + au.r + az.té + as + as” + as.f 
Nº de somas = 5 

Nº de produtos = 15 

Para Pn(r): 

Nº de somas = n 


Nº de produtos =1+2+3+4+5+..+n 


Ps(r) = ao + r(aj + r(a> + r(as+ r( as + às)))) 
Nº de produtos = 5 


E 
2 


Esquema de Ruffini-Horner ou Divisão Sintética 


ha 


- Pn(x) = Qua(x)(x-n)+R 


para x =: 
Pn(r) = Qua(r)(r-)+R 
.Pil)=R 


- Pa) = Qua 69(x-n)+ Qua(x) 
Pn(r) = Qua ((r-r)+ Qua(r) 
Pr(r) = Qna(r) 


DA PWN 


bs bs ba bz bi 


Cs.r Ca.r 


Togr=Pn(r) 


Exemplo: 
Ps(x) = 8 - 2x + 4x2 - 7x) + 5x" + x 


1 5 -7 - -2 8 
2 2 14 14 36 68 

1 7 7 18 34 76 
2 2 18 50 136 

1 9 25 68 170 
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Ps(2) = 76 
Ps'(2) = 170 
Fórmula Recursiva: 


bn = dn 

b; = Dj r+a;, j =n-1,n-2,n-3,....,4,3,2,1,0 
P(r) = bo 

cn= bn 

Cj=Cyl+ b; ' j =n-1,n-2,n-3,...,,4,3,2,1 
P(r)=c; 


Localização de zeros 
Seja a função f(x) e seja um x suficientemente pequeno. 
- 00 0) +00 
O problema de localização é encontrar um intervalo | = [a,b], tal que f(a).f(b) < O. 
Raízes positivas: 
f(O), f(Ax), (24x), ..., T(KAX), ... 
Testar se f(iAx).f((i+1) Ax) < 0 > i =[iAx, (+41) Ax]. 


Pode-se saber o número de raízes reais de um polinômio: 
Regra de Descartes 


Exemplos: 

Pa(x) = 1+ 3x — 5x2 +4xº + 8xº 

Nº de raízes reais positivas: 1 + 1 =2 
Pa(x) = 1- 3x — 5x2 -4xº + 8x” 

Nº de raízes reais negativas: 1 + 1 =2 


4.6 Atividades 


1. Verificar que o método de Ruffini Horner tem complexidade linear para encontrar 
o valor numérico de um Polinômio de grau n 


2. A equação x? -10cosx = O tem duas soluções: +1,3793646. Utilize o método de 


Newton para encontrar as soluções aproximadas, com precisão de 10º, usando 
valores iniciais xo iguais a -100, -50, -25, 25, 50, 100 
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4.6 Referências 
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Interpolação Polinomial 


META 


Resolver o problema: dada a função f(x), contínua ou 
um conjunto de pontos, aproximá-la por um polinômio 
de grau n. 


OBJETIVOS 


Estudar os principais algoritmos de construção destes 
polinômios. 
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5.1 Introdução 


Seja uma função f(x), contínua, uma das idéias mais antigas em cálculo 
numérico é aproximar esta função por um polinômio. 

Um Polinômio é fácil de manipular, encontrar suas derivadas, integrais e 
suas raízes com relativa facilidade. 

O teorema de Weierstrass afirma que “Toda função contínua pode ser 
arbitrariamente aproximada por um polinômio.” 

Os métodos a serem estudados como uma aproximação para uma função 
f(x) poderão ser aplicados quando, não conhecemos a função, apenas sabemos os 
pontos Xo, X1,X2,X3,............. Xn. Situação muito comum na prática quando se trabalha 
com dados experimentais. 

Se os pontos do parágrafo anterior forem distintos, determina-se um 
polinômio Pn(x) de grau no máximo n, tal que 

Seja o conjunto de pontos (x1, y1), (X2, Y2), (X3; Y3),..., (Xn; Yn) 

O problema de interpolação é encontrar um x para um y e [Xo; Xn]. 


5.1 Interpolação Linear 


Seja o par de pontos (xo, Yo)(X1, Y1), a equação da reta que passa pelos 
pontos é: 


Yo = dg + MXo 


= ao + ax, tal ve [ o 
y q Y1 = do + 0X 


h A Yo sd 
—1Lym “Dio ym 
(ani “41 Xo M 41 Xo 

1 X4 Fi a 


5.2 Interpolação Quadrática 


Para 3 pontos (xo, Yo), (X1, Y1); (X2, Y2) não-colineares. 
Seja P>(x) polinômio de grau 2: 


Po(x) = do + ax + ax? 

Yo = P>(Xo) = do + diXo + axo? 
y1 = Po(x:) = do + ajxi + azxi 
y2 = Po(xo) = do + axo + axo? 


2 
1 xo Xo Ag Yo 
1 x x ||jU|=|Y 
1 x, a a, y2 
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5.3 Interpolação para um polinômio de grau n 


Para n+1 pontos (xo, Yo), (X1,; Y1); (X2, Y2),... ;(Xn; Yn) não-colineares. 
Seja Pn(x) polinômio de grau n: 


Pn(x) = do + ax + axX+ ... + ax” 
Yo = Pn(Xo) = do + diXo + doXo + + AnXo” 
P 


n(X1) = do + diX1 + E Ro an anxi” 


Yn = Pn(Xn) = do + d4Xn + E: o quai? E EE qu 


1 x92.. x9] [Oo Yo 
de e fx = Vi 
RR RR Yn 


O problema de encontrar o polinômio que passe pelos pontos dados 
equivale a resolver o sistema de n+1 equações com n+1 incógnitas. Estudaremos 
métodos que resolvem esta situação em forma implícita. 


5.4 Método de Lagrange 


Este método é construtivo que engenhosamente pensó Lagrange, e que 
tentaremos reproduzir supondo quatro pontos dados ... 
Sejam os pontos . (Xo, Yo), (X1, Y1); (X2; Y2), (X3, Y3). 


Etapa 1 
Seja o polinômio de grau 3 construído da forma seguinte: 


Pa(x) = yoLo(x) + y1L1(X) + y2Lo(x) + Yy3La(x) 


Este polinômio será de grau 3 só se Lo(x), Li(x), Lo(x) e La(x) forem 
polinômios de grau 3, estes polinômios cnamaremos de Polinômios de Lagrange. 


Etapa 2 

O polinômio deve passar pelos pontos dados, ou seja, Ps(Xo) = Yo, Ps(X1) = Y1, 
P3(x2) =Y26 Pa(xa) =Ys. 

Para isto acontecer: 


Ps(xXo) =Yo, Lo(xXo) = 1, L1(Xo) = 0, L>(Xo) =0e La(Xo) =0 
Ps(x1) =Y4, Lo(x1) = 0, L1(x1) = 1, Lo(x1) =0e La(x1) =0 
Ps(xo) =Ya, Lo(x2) = 0, Li(x2) = 0, Lo(x2) =1e La(x2) =0 
Ps(xa) =Ys, Lo(x3) = 0, L1(xX3) = 0, Lo(xX3) =0e La(x3) =1 


1 izj 
ly) = to aj du 
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Etapa 3 
Os Li(x) são polinômios de grau 3. 


e plane a) 


EDGE) 
Li(x) = A 
Lo(x) = pm 
a (x — xo)(x — xx — x2) 


(x3 — Xo)(X3 — x9)(x3 — x2) 


Podemos verificar que depois destas três etapas o polinômio pode ser 
encontrado sem ter que resolver o sistema 4x4. 


5.5 Fórmula Geral 


Em geral para (xi,yi) =0,1,2,3,4,5, ...,n 
Teríamos que resolver um sistema de n+1,n+1 
A fórmula geral é dada pelas equações seguintes: 


Bi) =D y-Li08) 
i=0 
E (= Xj) 


js (Xo — Xj) 


ij 


Li(x) = 


Exemplo: 
Determinar o polinômio que passe por: (0,3) (1,5) (3,7) (4,9) e estimar o valor 
de y quando x=2 


Solução: 
Pa(x) = 3Lo(x) + 5Li(x) + 7Lo(x) + 9Lo(x) 
ger CC cs co 


(0-1)(0O-3)(0—- 49) - 12 
 GK-0k-3(x—- 4)  (x—0)(x—3)(x — 4) 
LO =q-)d- Dao” 6 
ee le Dr (ur 0-4) 
LO =3-)6-D0-9"” 6 
 G-0OG-D(x-3  (x-0)(x— Dx —3) 
Lt) = aa 4-3 12 
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3(x— Dx — 3x — 4) 5(x—0)(x — 3)(x — 4) 

Pa(x) E E o Vo BSS 
T(x— 0)(x— Dx — 4) 9x — 0)(x — Dlx — 3) 
POR O q 


3(D( DZ) ga (DD TOM) E, 90)(1)(—1) 


P;(2) = 12 


6 
1 10 14 3 24 É lê 
2'9' 92 2" 32 27 = 


5.6 Algoritmo de Lagrange 
O algoritmo é a implementação lógica das duas fórmulas dadas em 5.5 


P1. Fornecer os valores de (xi,yi) i=0,1,2,3,4,5............ hn, e o valor a interpolar x*, 
verificar se xi /xj 
P2. soma — O 
P3. Parai=0 atén 
prod — 1 
para j=0 até n 
Seif) 
prod — prod*(x*-xj)/(xi-xj) 
fim se 
fim para 
soma — prod*yi 
fim para 
P4. Mostrar soma // é o valor interpolado 


5.7 Programa de Lagrange no Scilab 


O programa foi feito para mostrar o polinômio e calcular após o valor 
interpolado: 


// Lagrange 

n=input('numero de pontos :'); 

[xJ=input('Digite os valores de x(i),i=1 n entre colchetes:'); 
[y=input('Digite os valores de y(i),i=1 n entre clochetes:'); 


for i=1:n 
for j=1:n 
ifi<>j 
if x(i)==x()) 
abort 
end 
end 
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end 
end 


xb=poly(0,"x"); 


yb=0; 
for i= t:n 
p=1; 
for j=1:n 
ifi<> j then 
p=p*(xb-x()V(x(i)-x(1)); 
end 
end 
Yb=yb+p“y(i); 
end 
yb 
xp=input('valor a interpolar :'); 
horner(yb,xp) 


Outros métodos que veremos a seguir baseiam-se no fato dos pontos 
estarem igualmente espaçados, isto é xX4-X0=X2-X1=Xg-X2= ........ =Xn Xni=h , € 
necessitamos definir um operador que facilite a notação das fórmulas que encontrarão 
os polinômios. 


5.8 Diferencias Finitas 


Seja o conjunto de pontos (xo, Yo), (x1, Y1), (X2 Y2),..., (Xn, Yn), tal que 
X=Xo+i.h, parai = 0, 1,2,...;n. 
Definimos o operador diferença A incremento h, como: 


ár) = a -2x+D=[K+h3-2g+hD+D]-[x3-2x+1] 


=xº+3x2h+3xh2+hº-2x-2h+1-xº)+2x—1 
=3x2h+3xh? +hº — 2h 


5.9 Propiedades do operador 


PI: IG) + 969) = Árto + Lg(0) 
pa: 4 = O, cconstante 

4 4 
P3: Left) = ci fm) 
Patos so 


ARA EA + 
Po: A h fo) =, f(x), menezZ 
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n 
P6:4 P(x)=c, cconstante 


AN+1 
P7: A Elio 


5.10 Tabela de diferencias 


2 AP yo —> Termo geral 


Exemplo: 
(0, 1) (1,2) (2,9) (3, 28) 


Solução: 


Termo Geral da Tabela 


Ali PA is 
h Yo» h Y; 
47, =Y — Yo 
42y0 = 4(4y9) = 4(y, — Yo) = 474 — 4y0 =(Y>—Y) — (Wi — Yo) 
=) —2% + 


4 yo = A(4?y0) = A(y, — 271 + Yo) = 472 — 247, — 40 
=(y3-y)-2(y —y) +04 —Yo)=Y3— 37% +34 — 
4ºyo = Ya — 43 + 67, — 494 +90 
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n 


yu =5 (ai! 


i=0 


5.11 Atividades 


1. Determine o tamanho do h = x; - Xi 


para a construção da tabela de f(x) = e* 
em [0,1] para que o erro de truncamento na interpolação linear seja menor que 
0.005 
2. Dado f(x) = sen x, f(0.1) = 0.09983 ; f(0.2) = 0.19867 Determine o valor f(0.16) e 
calcule o erro de truncamento 


COL Fr) = max t) FO)I) 


x 


3. Para f(x) = 5“ obtenha f(0.3) e o erro de truncamento se f(0.5) = 2.23608 


5.11 Referências 
CUNHA, Cristina. Métodos Numéricos. 2º Ed. Campinas SP: Editora da UNICAMP, 
2008. ISBN: 85-268-0636-X , CDD — 620.00151 


BURDEN, L. Richard, J. Douglas Faires Análise Numérica SP: Editora Pioneira 
Thomson Learning, 2008. ISBN 85-221-0297-X CDD - 515 
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Interpolação Polinomial 


META 


Resolver o problema de interpolação para pontos 
igualmente espaçados, gerando um polinômio de grau n. 


OBJETIVOS 


Estudar os algoritmos de Newton para a construção 
destes polinômios. 


linoOufs.br 


CÁLCULO NUMÉRICO Aula 6 


6.1 Introdução 


Os métodos seguintes usam as diferenças finitas na sua estrutura. Portanto, 
os pontos devem estar igualmente espaçados. 


6.2 O método de Newton para interpolação 


Seja o conjunto de pontos (xo, Yo), (X1, Y1), (X2, Y2),..., (Xn; Yn) igualmente 
espaçados: 
X1 = Xo + i.h, i =0,1,º, PER (| 


O polinômio de interpolação de Newton de grau n que passa pelos pontos 
dados é: 


Pn(x) = ao + ar(X - Xo) + az(x — xo)(x — x1) + as(x — xo)(x — x4)(x —x2) +... + 
an(x — Xo)(X — x4)(x — X2)... (X — Xn.4) 


X1) =Y1= ão+ as(X-Xo) 


Pn(Xi 

Pn(Xo) = Yo = do 

Pn( 

Pn(Xo) =Y>2=ao+ ar(x = Xo) + az(x = Xo)(x — X1) 


Ph(Xn) ii = ão + ai(X- Xo) + ax — Xo)(x — x1) +... + an(X — Xo)(X — X4)(X — xo)... 
(X— Xn1) 


Yi=Yo+as. h-> q = 
Y2o=Yo+ a,.2h + a>.2h.h 


= Quvo) 1—)o) — Yo 
h h 


— -2Y1 ++ A2 
p=yos e) oh +a>.2hº > q, = atvotaa — 1% 


2h? 2h? 
= Aéyo 
às gh 
= A'yo 
“Sun 


ata O usos = 


31h 


Ce) (6 = lr do) (E 21) 


Ayo Nº ya 
P. = +ou) + a (a 


n 


A”y 
teto 
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6.3 Notação fatorial decrescente 


(x — xo)(x — x4) = (x — xo)? 


Significa que tem-se dois fatores de base (x - xo) decrementando o outro 
fator em h. 


(x — xo)(x — Xo- ) = (x — xo)(x — x1) 


Em geral: 


xP =x(x-h)(x-2h)(x=3h)..(x=(n-1)h) 


6.4 Primeira Fórmula de Newton 


A primeira fórmula de Newton pode ser escrita em forma compacta usando a 
notação fatorial geral decrescente. 


Primeira Fórmula de Newton ou Newton 
Progressiva 


Exemplo: 
x y 
0 1 
1 ps 
2 9 
3 28 
h=1 


1 6 6 
P(x) = (XX) + smp (= Xo XX) + as (X— Xo MEX — Xi) (X — Xo) 
P(xy)=1+x+3x(x- D+ x(x- (x-2)=1+x+3X -3x+xXº 3x] +2x= 
=x) +1 


6.5 Segunda fórmula de Newton 
Seja o conjunto de pontos igualmente espaçados: 
(Xo:Y0) (XY) (Xos Yo) «+. (XmYn) 
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tais que: 


x=X+ih , i=1,2,3.,n 


Pa(x) = ao +a(x—x,) + as(x— x, )(X— X5.1) + as (X— X,)(X— Xn.4)(X— Xn.9) + 
+era (XxX, )(Xx— x, 4)(X— Xn.5)...(X— X;) 


P(x,) = y, vi 

[P(x o 

Ni n4 = +A(X,4—Xn) 

| 

IP(X,) = o. a = AG = Ma) 
nu Xn)(Xo — Xn- q «(Xo —X5) 

dq =Y 


Yn- =a Yn e Yn- AY n- 
Yna = do +a,(-h) —S a, = e 1 0 1 1 


h E h “h 
Yn2=ao +a;(-h)+a,(-2h)(-h) 
Ryr 
de= apê 
A Yi 
“= "ah 
a AY na es 
P(x) Yn + 1h (x— n) + 21h? (x—x,)(Xx— X,)+ + 
A"y 
+ jo XXX — 1) (Xx-x,) 


6.6 Notação fatorial crescente 
(xx, x —X04) = (x—x5)2 = (xx )x— x, +h) 


Em geral: 


x! = (x + hXx + 2hXx 43h). (x + (n— Dh) 
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6.7 Método de Aitken 


Este é um outro método para encontrar um polinômio que passe pelos 
pontos dados. E estes podem estar desigualmente espaçados. 

Seja o conjunto de pontos (Xg;Yo) (X,Y1) (Xos Yo) -.: (Xn5Yn), não 
necessariamente igualmente espaçados. 

O polinômio linear interpolante para o par de pontos (x9,Y9), (X,.Y1) é: 


Xx=x Xx=x 
PlO)=yo4 a E (Fórmula de Lagrange) 
0 1 1 
PAN Rae 
P(O)=Yo tota Xo) (Fórmula de Newton) 
1 Xx=x 
P()= o Yo 
X,—Xo XX Yi 


e Pos(x) 
Xi yí X-X1 Pora(x) 
> Prolx) o 
X2 y2 X-X2 2 Pogos..n(X) 
Piza(x) 
Ph n(x) 
Di n(x) 
Xn Yn X-Xn 


Po(x) é um polinômio de grau 1 que passa pelos pontos 
(X9,Y0) + (XY). 

Pi(x) é um polinômio de grau 1 que passa pelos pontos 
(XY) (Xp Yo). 
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P pie *=%Xo Yo PolXo)=Yo 
É Apos | Po(X)=Y4 
Pol) 1 XX Y Po(x,))=yY; 
so X)—X | XX Yo Po(X2)=Yo> 


O polinômio Po1>(x) definido como: 


1 
X> —Xo 


x—-Xo Po(x) 


Porz(X) = XxX, Po(x) 


é um polinômio de grau 2 e que passa pelos pontos (xo; Yo) (X1, Y1) (X2; Y2). 


1 
Po lXo) = [Oo — X0).Pio (o) — (Xg — X,). Po (xo )] = Flo) = Do 
AX — A 
1 
Fon (Mi) = [Ox — x9).Po (00) — (4, — x). PylxD)] = 
A — 
1 2h 
— h.y, —(—h = = 
es [h.y,—(-h)y,] hM 
1 
Pon (X) = [x — X9).Pp (x,) — (x, — x). By (x, )] = 
A —Xg 
1 
= (x, —X9).), =)» 
A —Xg 
1 x—Xo  Porz.na(X) 
P, Xx)= 
or2s.n(X) Xn—Xo0| XX, Prog na(X) 


6.8 Interpolação Inversa 


Seja o conjunto de pontos (x5.,Y0) (XY) (X5,Y5) «.. (X,Yn). 
Dado um y determinar o x. 
O polinômio deve passar pelos pontos (y,,X,) (Y5,X5) ... (Yo Xn) » isto é, 


Solução por Lagrange: 


linoOufs.br 


CÁLCULO NUMÉRICO Aula 6 


6.9 Atividades 


1. Determine : 


a) A: x b) AN Eua 
1 h 

Cc) AN sen x d) AN £ (x+h) 
1 


1 


e) AN x! £) Á «x + n) Em 


2. Encontre a formula geral para um elemento da tabela de diferenças finitas. 
3. Determine log 4.5 da tabela a seguir pelo método de Aitken 


log x 0.60206 0.62325 0.64345 0.66276 0.68124 


6.10 Referências 


CUNHA, Cristina. Métodos Numéricos. 2º Ed. Campinas SP: Editora da UNICAMP, 
2008. ISBN: 85-268-0636-X, CDD — 620.00151 


BURDEN, L. Richard, J. Douglas Faires Análise Numérica SP: Editora Pioneira 
Thomson Learning, 20083. ISBN 85-221-0297-X CDD - 515 
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Aproximação por Mínimos 
Quadrados 


META 


Resolver o problema de aproximação usando métodos 
de otimização. 


OBJETIVOS 


Estudar os algoritmos de Mínimos quadrados para 
diferentes tipos de funções. 
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7.1 Introdução 


A aproximação por mínimos quadrados é um método de otimização. Dados 
um conjunto de pontos (x; yi) =0,1,2,3,4.......... sn, a priori é definida uma função que 
tende a aproximar os pontos dados. Pode ser um polinômio, uma função logarítmica, 
exponencial ou trigonométrica. Escolhe-se uma métrica que meça os pontos dados à 
função, e escolhemos os parâmetros da melhor função que se ajusta aos pontos. 


7.2 Método de Mínimos Quadrados 


Seja o conjunto de pontos (x9,Y0) (X,Y4) (X5.Y5) -.. (X,:Yn)- 
Seja uma função f(x) que ajustará os pontos dados, através de uma métrica d,. 


Exemplos de métricas: 


d, =f(x,)-—); 
Id HfG)->,| 
(d;)? =(f(1)-9) 


min 5a, 
min bo CA 
1=0 
Mínimos Quadrados 
Seja f(x) = P(x): 


Pl)=a+a,.x 
d, = PR(x,)— 9; 


à Oba 


n n 
min 3 d?=min o (ay +a x, — 7) 
i=0 i=0 


G(aç.a,) => (ay + ax, —);) 


i=0 


Kaga) 6 
Po 

Aaga) 6 
du, 
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n 


Ha, +a,.x, -y;)(D=0 
i=0 


1 


n 


2 ay +a,.x, -y;Xx,)=0 


i=0 
ba A, +aD, X; -5. y,=0 
a X +aD, da -5. 194, =0 
(n+Da, +) x = y; 
>, A +a,;>) me = Rad 
ba ; » A; n+1 Z Y; 
E > Xy :» a ba X; > Ae); 
“o I(n+] mM (n+D) Dx; 
2% 2a >, 4 DE 


Se o polinômio for de grau 2: 


Blw)=atax+a,x 
d, =P(x,))-); 
n 
E 2 2 
min z (a; PAX Fa 9) 
1i=0 
2 


n 
2 
Gaçsa,a,)=> (agtax, tax) —y;) 
1=0 


DS Aa +rax, tax) -y)D=0 


> Ma, +a,x, O PR —y)(x,))=0 


> Xaç+ax, tax -y)x) =0 


i=0 


+; ZA >. la, |= 3x); PAi Matriz Simétrica 


n+1 bo Ed j bão 
sa pia Da Ma 


A, 
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Em geral: 


n+10 3a, bx Sax la 2 
3x, sa Dx) E jul. as 


n n+l n+2 n+n n 
Dx Dx Dao Dx] a) Dx, 


7.3 Atividades 


1. O volume de álcool anídrico em função da temperatura esta dado pela tabela 
abaixo: 


Temp (Graus C) 13.9 43.0 67.8 89.0 99.2 
Volume(cm3) 1.04 1.12 1.19 1.24 1.27 


Fazer um ajuste para v(t) = 1 + bt + ct? 
Construir a tabela v = v(t) para t = 20(5)40 


2. Dada a tabela 
x 1.0 1.05 11 1.15 1.2 1.25 1.3 1.35 
y 1.0 1.01 1.02 1.04 1.05 1.06 1.065 1.08 


Estimar f(1.22) por regresão linear. 


7.3 Referências 


CUNHA, Cristina. Métodos Numéricos. 2º Ed. Campinas SP: Editora da UNICAMP, 
2008. ISBN: 85-268-0636-X, CDD — 620.00151 


BURDEN, L. Richard, J. Douglas Faires Análise Numérica SP: Editora Pioneira 
Thomson Learning, 20083. ISBN 85-221-0297-X CDD - 515 
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Integração Numérica 


META 


Resolver uma integral usando aproximação polinomial. 


OBJETIVOS 


Estudar os algoritmos que resolvem em forma 
aproximativa a integral de uma função e estimar o seu 
erro. 
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8.1 Introdução 


Os métodos de aproximação polinomial são usados para integrar 
numericamente uma função y=f(x) num intervalo dado [a,b] ou mesmo um conjunto de 
pontos (xi, f(y;)) =0,1,2,3,4,......... hn. 

Casos em que a função é difícil integral ou não tem solução analítica, um 
polinômio sempre é de integração imediata. 


f(x) 


A área fechada em vermelho representa a integral definida do polinômio, e a 
línha em azul é a função. 


8.2 Integração Numérica 


Seja a integral definida da função f(x): 
b 
|roras 


A integração numérica é utilizada quando não conhecemos a função, e sim 
pontos dela, ou a função não é uma função integrável analiticamente. 


Solução: 
b b 
[food = [Pax 
Aproximação Linear: 
b 
Ia, +a s)dx=[a,+ax/2) =a,b+(a,b'/2D)-aça-(a a? (2) 
Dados dois pontos de f(x): (X9.Yo) € (x,,Y 4). 


linoOufs.br 


CÁLCULO NUMÉRICO Aula 8 


P(X)=Yo+ AFo fg X,)) (Newton Progressivo) 


h 
onde : 
h=(x,—Xo) 
Yo=Y1—Yo 
Edo é 
D=Yo 
Jrcoar- fio cê Xo)ldx = Yo “= 
Ay Áyo X 
= pat [2-X%,.X9)— Y9Xp — E Do w)= 


= Yo(X, — 19) + (o [ar do Aa “12)= 


Ayo 1 


— Jo: 


h 
Jrcode= [o +21] 


A=[yo + v15 ; A= Área do Trapézio 


Para (X9,Y0) (X,Y4) (Xp Yo) : 


Trogax = que 


Xo Xo 


Jree- Jon ema A DR eras 


h 


X, 


Seja t= 
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Para x=X, > E CE e 
o XX 0 
Para x=x, > t= R 0 
(x-x,)2 =(x-x,)(x— x, +h) = 
=(x—xo)(x— (x, —h)) = 
=(X—X,)(X— x) 
(x=x0)2 (x-x5) (x-x5) 
dE RCE” =t(t+1) 


dt= ds > dx=hdt 


0 2 2 3 2 
Ao Ayit A Yo t t 0 
+ Ayit+ t(t + Dlhdt = hly,t + + + = 
Jo» it + SG +D] De 
A yo [-8 4 
=h[2y, -2Ay, — disp li= 
[27, 4 2 [5 5) 
Ay, 2 
=h[2y, —2Ay, — O. = 
[2% y, 2 3! 


A? 
=h[2y, — 24y, — o 
l 2 
o, —6Ay, A YolJ= 
1 
= om —6y, +6), - NA yol= 


1 1 
= oo +) —2% + Yol = h6», +) 2%) + Yo] 


8.3 Fórmula Geral (Newton - Cotes) 


Seja o conjunto de pontos (X9,Y9) (XY) (X5,Y5) -.. (X,5 Yh). 


lrogdx = PP o)dx 


P(x)=> y-L(X) (Fórmula de Lagrange) 
= 0 


Sejam os pontos igualmente espaçados: 
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L()= Mem Ada las) 
O us sa =) sms Gs (u=m) 


Se x=Xp+ih: 


Lis G=mÃe-a reader) ESA) 
j ih(i—Dh(i—2).h.....h(—Ah)(-2h)....—(n—iDh 


X—X 
Seja S=——, então: 
ORA A 
h h h 
TA LEA, 
h h 


Liss S(S —D(S — DS — 3). S —(— DS — (+). dS —n) 
iG -1.6-2)..1(D(C2)....f—1) 

ng CD. 
Mn -i)MS —i) 


ds =x — dx=h.dS 


E nn n n dd (—1) n=i 
h |P(S)dS= h “Lo(S)dS=h fo SE. 
A ps e O nO 


Sen=1 -»> Trapezoidal 
n=2 > Simpson 
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n 1 SO nr 1 gi 1 e (E 
h5 y |——— ds =h d —— 
2, RASDNSTO i Do ma cos=O) di A SDISID 

ES SS 


ds]= 


1 1 
=hyo Jos —I-Dds+y, [Sds]= 
0 


1 


Ss? 2 


1 
1 S 
= hl- yo 08h + = 
2 ho 


l 1 
= hyol; Mt nlçll= 


l ) 
Hd Ei Dio + 


8.4 Método de Romberg para Integrações Numéricas 


Seja o conjunto de pontos (xp, yo XX, Y XX, VD... (X,,Y,) 


Irtgdx 


A idéia de Romberg é repetir fórmulas que implicitamente geram polinômios 
de interpolação de grau n. 


b a+hy 


[rooax — [rodax + [rá 


a a a+h 


a+h, a+2h, a+3h, 
Í fO)dx= | fO)dx + Í fO)dx + Í fO)dx + I fO)dx 
a+h, a+2h, a+3h, 


2" a+(+Dh, 


frovas= 5 [rosas 


i=0 a+ih, 


As integrais são aproximadas pela Trapezoidal: 
: h 
n=0 > [rod h=b-a 5 T E te O] 


n=1 5 T = fa) + Nash) Hy(a sh) + fOI= 
É = fa) +2f(a+h)+ 1] e 
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n=2.5 T = ELA) + fla + ho ELI (a ho) + Ha + 2h) 
+EL/ a +2h,) + Fla + 3h, JH ELA (A+ 3h) + FCO), h, + 


hy aê 
n=k > Tess ie NG) O] 


Para entender a idéia de Romberg é necessário saber o erro na fórmula 
Trapezoidal. 


[rodas a ALf(a)+ O) h=b-a 
ffeodr=V. e Trt)+fO=V, 


es =Vo-V.= [Gde SL ta) + f(b) 


a+h 


et)= [Hdr ft) + Fa +] 
cn= flash) -SL Fla += (a) + fla +h) 
E pd oh n A ' ad ' 
e'r= far far farh)-— Sah) 


Er= Sa +h) 
1 =max(f(a +) 


e" (h)< 


h " 
5! (9) 
e”, (h) < 1") 
" h ” 
forr dh< [5 109 dh 
e (h)< "0 dê 


Seh=0: 
e(0)<0+c > c=0 


e(0)= f(a+0)-S-+Uf(a) + fla +0)]=0 
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1 ho "W 
Jean [-— p"eeydh 
hº 
elh) < ET (S)+e 
nº 
elh) < RÁ (5) 
O erro é calculado somente se conhecer a função f(x). 
' h 
|rodde=T, =1100)+ 10] 
h=b-a 


[redr=m =Lyta+2/0a+h)+ 10] 


KR 
Den o) 


ho " ho ” 
ut E ua! (55) 
3) f(S)=maxtf (o), FS) 


2h, 
ey =— E 116) 
ME ata DO ai 
Sjh= então £y, = õ f'(o) (1) 
8h, " 
6) e, =— 12 (6) (3) 
NDe(4)e (6) 
2h; 
Er =— E (e) 
4.2hº 
En=-q2 1 (O 
Ser =4er 
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Er, = fi (o)dx —T | Er Í fO)dx—T, 
[roa-t, =4[fcar-T]) 


4T,-Ty =3 | f)dx 


frear = Fido 


47, =4 Ea) +2/(a+h)+ FO 


2h 
E Lita) + 1(b=DLyta) + f(b] 


47, To = (HM) +2 Fa +h) + FD ha) + FD) = 
=h(f (a) +4f(a +) + f(b) 


4T, —T, 


h 
3 E EA 


Fórmula de Simpson 


8.5Tabela de Romberg 
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4To = Tê 
T = 1 a 0 
4TE=T 
E = E 1 
AT = Tº 
1 [E i 
=" 
16T)-=T; 
2 1 0 
tj= 15 
Sn .=d" 
2 1 i 
ti 15 
4º T? o T? 
E 
Em geral: 
rara 
jo [ia i 
f="4 
8.6 Atividades 


1. Determinar fórmulas para integrar 


[e fl)dx e po f(x)dx usando aproximações de f(x) por polinômios de 


interpolação Newton progressivo e regressivo. 
2. Calcular a integral í e “dx utilizando a fórmula trapezoidal para n=10 e estimar o 
erro. 


3. Calcular as seguintes integrais pela fórmula trapezoidal e Simpsom com erro 
menor que 0.01. Determine o h que faz o erro menor que 0.01 


a) [axid+x) b) frinadr 
o) o) 

(o) [ex d) [cosx / xdx 

4. Encontre a fórmula geral para a regra trapezoidal n intervalos igualmente 
espaçados 


5. Seja o intervalo ho = b-a, hy = ho/2,...... hm = hn4/2 Encontre as trapezoidais TT, 
sas ii 


6. Determine a fórmula de T' em função de T' 
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8.7 Referências 
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Solução de Sistemas Lineares 


META 


Resolver o problema de equações lineares de qualquer 
tamanho. 


OBJETIVOS 


Estudar os diversos algoritmos, analíticos e 
aproximativos e sua implementação no computador. 
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9.1 Introdução 


Muitos problemas de engenharia e pesquisa operacional são resolvidos 
usando a álgebra linear. Isto é, matematicamente são reduzidos estes problemas a 
um sistema de equações lineares. Por exemplo: Cálculo da tensão em estruturas da 
construção civil, solução de equações diferenciais parciais, determinar o potencial em 
redes elétricas, problemas de otimização, etc. 

Quando o sistema é de grande porte, devemos ter cuidado de preservar ao 
máximo a melhor exatidão e precisão. 


9.2 Solução de Sistemas Lineares 


Seja o sistema: 


> 
x 
II 
o 


onde: 


A=(a;;)am + X = (x,),b=(b;) 


i=0,1,2,..,n 
p=" 0, 1,2, 
OU 
A dy Gn || b, 
dy dy dy don ||Ão b, 
dy dy das An || X3 |=|D; 
di E Cx An | | Xn bn 
OU 


AX, +a,,X, + AX +.. tax, =D, 

AX, + AooX, + AogX4 +... +AGnX, = Do 

AsXj + Aa9X; + Ag Xg +... +AGnX, = Ds 

AX; + AoXo + An5X +... +AnX, =D 
ou 


DS oa ea E 


/ 
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Um sistema linear nxn que admite uma única solução é chamado de 
determinado, se admite várias soluções é dito de indeterminado, e se não admite 
solução ele é impossível. 


9.3 Solução algébrica: 
Ax=b 
1) Se o determinante de |A|x 0, então existe inversa da matriz A, A”. 
2) Multiplicando a esquerda por A”: 


AT.Ax=A".b 
Ix=A.b 


3) x= Ab (Solução teórica) 


Na prática, se o sistema for de ordem n>5, há dificuldade de resolver em 
forma manual. 


9.4 Método de Eliminação Gaussiana 
Seja o sistema: 


dj), + djoX + AX + AX, = djs 
dpyX, + dooXy + dogXa + dogXy = dos 
Ag1Xy + dooXy + EosXs + AX, = dos 
AX, + daoXo + dgsXa + AX = gs 


O método consiste em transformar o sistema Ax=b em outro sistema 
equivalente Dx = f, tal que, D é uma matriz triangular superior. 

Para isto, utilizam-se as propriedades das equações: 

P1: Se multiplicamos por uma constante uma equação a equação não varia. 

P2: A soma de duas equações é linearmente dependente as equações 
somadas 

A transformação ocorre usando estas duas propriedades 

Ax=b ss Dx=f 


A matriz D resultante é triangular superior. 
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Exemplo de matriz triangular superior: 


212 3-1 
Dl =2 5 4 
D=|0 0 5 2 0 
00 0 1/1 
00 0 0 3 


9.5 Algoritmo de triangularização 
Passo 1: Sea, =0 então “Troca linha 1 por linha i, i = 2,3, 4º 
linha 1+<— linha 1 


linha 2 — - da) “linhal+ linha 2 


dA 


linha 3 [ du) linhal + linha3 


dA 


linha 4 — [ da “linhal+ linha 4 


dj 


Passo 2: Se a>>= O então “Troca linha 2 por linha i, i = 3, 4º 


linha 1 < linha 1 
linha2 < linha 2 


dao 


linha3 - a beimha 24 linha 3 


22 


do 


linha 4 e [8 beimha 24 linha 4 


22 


Passo 3: Se as3= O então “Troca linha 2 por linha i, i = 4º 


linha 1 < linha 1 
linha 2 < linha 2 
linha3 <— linha 3 


linha 4 € [Se ) tinha 3 inha 


33 
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Para encontrar o termo geral definimos três índices. 
Indice para o Passo: k = 1,2,3; 

Índice para a linha: i = k+1,....,4; 

Indice para a coluna: j = k.......... se 


Algoritmo: 


Parak = 1,2,3 
Se a = O então Rotina Troca 
Para i = k+1 até 4 
Para j=k até 5 


Fim 
Fim 
Fim 
Para qualquer N: 
k = 1 até N-1 
i=k+1 até N 
j=k até N+1 


X, = dus [au 

Xa = (das — 24X) / Ag 

X, = (Aos — do9Xa — AoaX4) / do» 

X,= (aj — AX, — AX — AX) / A 


Termo geral: 


Xn= anns1/ Ann 


N 
Ae (a, NH Zax,)/a; 


r=|+ 


Jeni squesraradas 3,21 


Exemplo: 
| 3X, +X =—2 
X+X%+X=3 
X—-2X,+X,=0 
Solução: 
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Passo 1: k=1 
O =3. 1 =2 111838 
1 1 1 93 Linha 2 «> Linha 1 0 -9 1 -2 
1.210 1210 
11183 11138 
0-3 1-2 Linha 3 < (Linha 1)(-1)+Linha 3 03 142 
1 210 Oo3S0Sl4 
Passo 2: k=2 
11138 11 1/3 
03 1-2 Linha 3 < (Linha 2)(-1)+Linha 3 03 49 
Oo30S34 Oo 0 +41 + 
PaiarçaÉ X,=1 
-SX+X=-20 > X, =1 
—X = —1 el 


9.6 Método de Gauss-Jordan 
Seja Ax=b. 


O método para a solução do sistema consiste em transformá-lo em outro 


sistema identidade Ix = b, usando as mesmas propriedades das equações aplicadas 
no método de triangularização. 


Para uma matriz 3x3: 
Passo 1: Sea,=0 => RotinadeTroca i=1,2,3 


linha 1 < linha 1/a,, 
linha2 <— (- a,;).linha 1 + linha 2 
linha3 < (- as;).linha 1 + linha 3 


Passo 2: Se a,=0 => Rotinade Troca i=2,3 


linha2 < linha 2/a,, 
linha 1 <— (- a,,).linha 2 + linha 1 
linha3 «<— (- as,).linha 2 + linha 3 
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Passo3: Se as=0 => Rotinade Troca i=3 


linha3 < linha 3/ ass 
linha 1 <— (-a;).linha 3 + linha 1 
linha2 < (- as).linha 3 + linha 2 


Passo >»k=1,2,3 
Linha >i=1,2,3 
Coluna > j=1,2,3,4 


Algoritmo: 


Parak = 1 até N 
Se a = O então Rotina Troca 
Parai=1 até N 
Se i=k então 
Para j=katé 5 
q; = d; [ A 
Fim 
Senão 
Para j = 1 até N+41 
dj = (ax). A +; 
Fim 
Fim 
Fim 
Fim 


Solução do sistema: 
RESET 2S a N 


Exemplo: 


pe =-2 
PETS 
X—-2X,+X, =0 


Solução: 
O =3 1 =2 RR TR TR 
1 1 1 93 Linha 1 «> Linha 2 0 = 1 2 
1210 Lo =2. 10 
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73 11 13 

0 — 1 2 Linha 3 < (Linha 1)(-1)+Linha3 0 -9 1 2 

1-2 10 O 30 Ss 

po tod 3 ; 1 3 

0 | 1 2 Linha 2 < Linha 2/(-3) 0 1 A / 93 -2 / 93 

O 303 O 3 0 —3 

1 1 1 3 Linha 1 <— Linha 2(-1)+Linha 1 1 0 4 / 3 7 / 3 
0 1 + / 93 -2 / 93 Linha 3 < Linha Pla) linha a 0 1 A / 93 —2 / 93 
O 3 0 —3 00 +41 — 
10 4/8 1/8 maicimasaaaima 10 01 

0 1 + / 93 -2 / 93 Linha 2 < Linha 3(1/3)+Linha 3 0 1 0 1 

0 0 + — 0 011 


XG=1,M=1,X=| 


9.7 Atividades 


1. Resolva o seguinte sistema de equações pelo método de eliminação gaussiana 
usando as funções do Scilab. 


X= = Zed w+ x+y + z2z=10 
5x - 4y + 32 =-12 2W+3Xx+Yy +92=831 
2x+ y+ z= 11 -w+ x-Dy+3z2=-2 


Sw + x+7y-22=18 


2x+6y- z=2 
5x- y+22=29 
-3x-4y + z=18 


2.- Resolver pelo método de eliminação gaussiana , método Gauss-Jordanl, o 
seguinte sistema tridiagonal ou matriz banda, usando as funções do Scilab. 


2x1 - x2 
-x1+2x2- x3 
- X2+2x3- x4 
- Xx3+2x4- x5 
- X4+2x5 - x6 
-x5 + 2x6 


| 
PM an — 


E VR Re GR 


ll 
pm 
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Solução de Sistemas Lineares 
(continuação) 


META 


Resolver o problema de um sistema linear, de qualquer 
tamanho. 


OBJETIVOS 


Estudar os algoritmos de fatoração LU e métodos 
iterativos. 
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10.1 Introdução 


Os métodos de fatoração são especialmente úteis quando se tem que a 
matriz A pode ser expressa em um produto de matrizes LU, onde L é uma matriz 
triangular inferior e U uma matriz triangular superior, definidas adiante. 

Se os valores iguais a 1 estão na diagonal L, o método é chamado de 
método de Doolittle e se os valores 1 estão na diagonal U , o método é chamado de 
método de Crout. 

Os métodos iterativos são aproximações sucessivas de vetores solução que 
tendem ao valor exato no limite. Requerem uma condição de convergência. 


10.2 Fatoração L.U. 


Seja o sistema Ax = b. 
O método consiste em transformar a matriz A em um produto de matrizes 
triangulares: 


A = LU 
A = (aj) rn 

[1 0 0... 0] 

n | 

| tor lo 0 | 
L | Ei. la lg . O | 

, Ino Ina N 

Ê U,> Uia A Ur | 

| 0 1 Usos Usn | 
U E o | 

lo O 0. 1 | 


Ax = b => LUx = b => fazendo Ux = z 


Lz = b Sistema triangular inferior resolvido em forma recursiva e: 
Ux = z outro sistema triangular resolvido recursivamente. 


a, d> an | [; 0 0 | [1 U,, vd] 

| || | | 

| dm dy don [=| lo; Io 0 || 0 1 Up | 

| NIE E e] 

ES ano a, | [1 a aaa || Lo 0 1 | 
A = L U 
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ly=as 
help =a Up =a/l,, 1,40 


“Multiplica cada linha de L com todas as colunas de U” 


10.3 Métodos Iterativos para a Solução de Ax = b 
Seja o sistema Ax = b. 


Passo O: 
Transformar o sistema Ax = b em outro sistema equivalente de forma : 


x=Cx+f 
Passo 1: 
Valores iniciais: 
X =f 
jo 
Passo 2: 
jej+1 
Ret sf 
Passo 83: 


Selx'-x"| <e vi 
então Solução aproximada x! 
Senão Volta ao Passo 2 


O algoritmo gera uma sequência (x) > x* como solução. 


10.4 Condições de Convergência 


Para que a sequência gerada, (X'), seja convergente, é necessário que 
|C| <1 (|IC| - norma da matriz). 
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Norma da matriz C 


IC] = mín |C PCS 
Cl, — norma da matriz linha 


Cl, = max Deal 
Cl. =max (Bleu 


Cl <1 


“Ou a soma de todos os elementos das linhas ou a soma de todos os 
elementos das colunas deve ser menor que 1”. 


Er 
Exemplo: A=1 0 1/2 0 
Ka 1/4 e 
Solução: 
Al =míntJAl],. A], 


Are) 
a (AQ) 


Seja o sistema Ax = b, sistema Diagonal dominante, então: 


AX, + dj2X> + dig Xa +. + Ankh = Ein 
dpyXy + dooXo + dogXg +. HdonX = don 


AX + dnoXo + dn3Xs +. +dnXn =— Annui 
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Exemplo: 
[Ba mtemão =] 
x, +4x, —x,=4 
X,—X,— 5X =5 
Solução: 
1 
X=0X+gX+g% + 


1 1 
Xo=—4% + 0x, + q +1 


1 1 
X, o 


Ixl lo 1/3 1sllx! [1/3 


ef 1/4 0 1/4 perl: 
x] [15 145 1 lx] li 


x=Cx+f 


EA 


= mín t|C 


cl; 


7? 


Q 
IH 
E 
pá 
p= 
REA 
[a 
+ 
O |a 
no— 
| 
O 
+ 
pim 
Na 
(TN 
njm 
| 
O O 
isa 
—— 
IH 
a SIN 


Q 

Il 

E 

Ss 
p=, 
WIN 
fm 
o | 
E 

A 

pa 


10.5 Método de Jacobi: 


[1/3] 
Passo 1: “| ç , €=0,1 


1 
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Passo 2: 


Passo 3: 
Regra de Parada 


<g 
Eis 
12 


18 
15 


1 
4.— 
3 


1 <g 


1 


<e 


10.6 Método de Gauss-Seidel: 


Passol: 
[1/3] 
X =| 1 
ES 
Passo 2: 
= 
1 1 
do — e — EEE = 
Ms ad a 1 
: 1 1 
Xp == "1+0x, + q d+1=1 
1 1 
= gas] 
10.7 Atividades 


1. Transformar as matrizes em fatores LU 


128 128 
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2 4 6 De dd É 

> 6 7 2 4 6 

2. Escrever o algoritmo, e as formulas gerais para o método do elemento maior que 
funciona igual ao método de Gauss-Jordan sendo que a escolha do elemento 


pivô é o maior elemento da coluna em valor absoluto, entre as linhas que não 
contém elementos pivôs escolhidos. Fazer trocas de linhas para arrumar o 


sistema. 
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